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ABSTRACT 


Radio  asti'onomy  star  tracking  requirements  are 
resolved  into  elevation  and  azimuth  cooi’dlnates  and 
related  to  servo  system  transfer  function  specifica¬ 
tions.  Time  functions  for  angulivr  position,  velocity, 
and  acceleration  are  derived  and  maximum  values  found. 

A  nijmerical  example  based  upon  the  latitude  of  Washing¬ 
ton,  D,  0,,  and  a  star  declination  such  that  at  transit 
the  angular  distance  from  the  zenith  is  10  degrees  or 
greater.  Is  used  In  the  calculation  of  the  azimuth  servo 
bandwidth  needed 'to  insure , dynamic  errors  appropriate  for 
automatic  tracking  with  a  60-foot  parabolic  reflector 
(tracking  error  less  than  l/4  degree ) .  A  bandwidth  of  a 
0,1  radian  per  second  is  shown  to  he  adequate.  A  block 
diagram  is  given  for  a  servo  system  suitable  for  the 
azimuth  qoordina'te. 
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ELEVATIOK-AZIMUTK  3SRV0  SYSTEM  SEECIFICATIONS 
FOR  STAR  TRACKING 


nmioDUCTiOK 


The  advent  of  radio  astrorion^  has  seen  employed  a  wide  variety  of 
antennas  and  means  for  beam  steering.  Among  these  is  the  parabolic  re¬ 
flector  mounted  on  a  pedestal  with  elevation  and  azimuth  coordinates 
of  motion.  Seryo  system  requirements  for  such  "EL-AZ"  mounts  become 
automatically  defined  upon  specification  of  the  tracking  problem'.  In 
the  study  reported  here,  the  tracking  problem,  a  time-domain  phenonemon, 
is  transformed  into  the  equivalent  frequency  domain  ftinctlon  so  that  the 
servo  systems  my  be  designed  ixsing  the  frequency  response  method.  In 
choosing  the  open-loop  ti-anafer  functions  ■  for  the  seryo  systems  several 
factors  are  considered,  naiEely:  (l)  close  fitting  of  the  excitation 
function  to  achieve  the  laaxtmum  exclusion  of  noise  simultaneous  with 
adequate  bandwdlth  for  transmission  of  frequency  components  in  the  input 
signal,  (2)  auilable  gain  aiid  phase  uiai'glns  'Lo  insui^e  closed- loop  servo 
system  stability,  and  (3)  the  proper  gain  to  insure  that  the  maximum 
dynamic  tracking  errors  almost  equal  but  do  not  exceed  the  specified 
limits. 

The  derivation  of  time  functions  for  elevation  and  azimuth  angular 
position,  velocity,  and  acceleration  for  star  tracking  is  presented  in 
the  initial  part  of  the  body  of  the  reiiort.  Gulllemin’s  (l)  method  of 
approximating  the  frequency  domain  transform  of  a  time  function  is  tnen 
applied  to  the  above-horizon  portion  of  the  azimuth  time  function.  A 
transfer  characteristic  Is  then  chosen  follo'vring  Graham's  (2)  procedure, 
in  which  the  steady-state  error  series  (3)  plays  a  central  i*ole.  The 
block  diagram  for  a  possible  mechanization  is  given.  The  trajectory 
equations  for  the  elevation  axis  are  xjrovided  so  that  a  similar  study 
for  the  elevation  coordinate  servo  nv>y  be  made  following  the  came  pro- 
eadui^e . 

CELESTIAL  AND  TERRESTRIAL  GOOKDiHATES 

Star  tracking  euagvilar  position,  velocity,  and  acceleration  quan¬ 
tities  resolved  into  eleyation  and  azimuth  coordimtes  are  desired  as 
a  function  of  time .  The  location  of  maximum  values  of  velocity  and 
acceleration  in  both  axes  is  also  needed.  In  the  derivations  and  in 
Fig,  1  a  modification  of  the  notation  used  by  Hosraer  is  employed. 
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Let 


s  =  interior  angle  of  aGtroRomical  triangle  at  zenith,  degrees 
T  ^  hour  angle  (=ajfc) 

CO  =  apparent  angolai'  velocity  of  the  stars  about  earth's  axis, 

in  radians  per  second  =  .2tf  =  7,272 *10" in 

24- 3600 

degrees  per  minute  =  l/4 
t  =  time,  seconds  or  degrees  of  rotation 
0  “  azimuth,  degrees  (0  =  -z) 
h  =  altitude  (elevation  angle),  degrees 
S  =  declination,  degrees 
4i  =  latitude,  degrees 


AZIMLfPH  ANGULAR  Pi'JSITION 

To  find  the  azimuthal  quantities,  start;  with  Eq.  (31)^  page  3^ 
of  reference  4, 


sm  z 


-  3m  T  Cos  S 


Cos  h 

Substitute  from  Eq.  (ll),  page  3?  of  refereiice  4, 

CoS  h  =  S  cos  ^  -  cos  S  sm  ^  cosT 

cos  z 

and  obtain 

sm  z  sin  T  cos  8 


sin  S  cos  4  “COS  S  sin  ^  cosT 


(1) 


(2) 


COS  z 


*l;dn  z  s 


_ _ sm  T 

Cos  ^  -  sin  A  e©sT 


(M 


Numerical  evaluation  of  Eq.,  (4)  gives  an  angle  z  which  is  iSO*^  from  the 
correct  angle.  This  occiu’red  through  the  use  of  the  sine  formula  with 
its  inherent  ambiguity.  Making  this  correction  and  the  changes  of 
variable  T  =  ojt  and,  z  =  -0,  we  have 


■tnn  6  s 


-  sin  {ji't 

A  -  B  cos  u>*t 


(5) 


where 

A  -  COS  d)  'trdvi  S 

8 

B  =  sm  4 

AZIMUTH  ANGUIAR  VELOCITY 

Azim’uth  angle  0  aa  given  by  Eq.  (5)  may  be  differentiated  to  obtain 
the  corresponding  angular  velocity.  Thus, 


see*©  is 

dt 


{A"8  eos  ■ngto 

(A«B  c®si8*t)* 


A  u)  eoA  0'i;  -t*  Bio 


(A-  B  cos  u)t)* 

rB-Acasu>tl  a 


0} 


sin*© 


(T) 


(8) 


Azimuthal  angular  velocity  is  given  by  Eq.  (8)  except  at  points  where 
sin  ojt  "  0.  For  these  points,  resturn  to  Eq.  (7)  and  write 
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(9) 


dt 


itsO 


_ to 

&in  ^  CAS  ^  tavi  3 


Figure  2  sticws  the  azimth  angular  velocity  as  calculated  for  t  =  0 
using  Eq.  (9)  and  a  latitude  4i  =  39°  (approximately  the  latitude  of 
Washington,  D.  G.)*  The  extreme  values  shown,  S  =  ^  -  10®  =  39°-  10*^  = 
49^^  and  29°,  are  hereafter  referred  to  as  Case  1  and  Case  2  respectively. 
Since  Case  2  exhibits  a  higher  absolute  value  of  azimuth  velocity  for  a 
given  approach  to  the  zenith,  detailed  calculations  are  made  for  Case  2, 
Returning  to  Bq,  (5),  the  constants  become 

Case  2 

4>  =  39'^ 

S  =  29° 

A  =  cos  39°  tan  29^^  =  O.77715 '0,55431  =  O.43078 
B  =  sin  39°  =  0.62932 

Figure  3  shows  Eq,  (5)  calculated  for  Case  2  over  a  24-hour  i)eriod. 
The  Indicated  points  corresponding  to  the  star  rising  and  setting  are 
from  calculations  shovm  later  in  the  report.  Using  Eq.  (9),  the  maxi- 
mxjm  azimuthal  angular  velocity  becomes 


Case  2 

de  l  .  7.273  '10^ 

dt  |t*0  0.62932  “0.4^070 


7.272  _ 

0.198S4 


0.3663 


radian  per  second 


Thus,  for  Case  2,  at  t  =  0  the  8.zimuth  angle  is  9  =  l80°  eind  the  angular 
velocity  is  0  ~  O.3663  rradian  per  second  =  1.259  degrees  per  minute. 

For  S  ^  i>  (a.s  in  Case  2),  the  azimuth  angle  increases  continuously  from 
star  rise  to  set.  For  3  >  <>  (as  in  Case  l),  the  azimuth  angular  rate 
is  initially  -pOiSitive  but  reverses  twice  before  the  star  sets. 


IMUTH  AVGLE 


Of  interest  in  servo  system  design  is  the  2.’elative  period,  during 
which  the  input  angular  velocity  is  less  than  the  liiininnuu  velocity  of 
the  servo  motor,  as  installed  in  the  mount.  Good  design  may  yield  a 
dynamic  range,  or  ratio  of  maximum  to  minimum  motor  speed,  as  high  as 
100  to  1;  thus  the  region  of  inputs  where  the  rate  is  less  than  1^ 
of  the  maximum  may  be  called  the  "near-zero  rate  region."  This  region 
is  related  to  star  tracking  in  the  azimuth  coordinate  from  a  station 
in  the  northern  hemisphere  as  follows: 

a.  Stars  with  south  declination  and  stars  with  north  declina¬ 
tion  less  than  the  latitude  of  the  tracking  station  exhibit 
cuiitinuous  azimuth  rates  from  east  through  south  to  west, 
with  minimum  values  at  star  rise  and  star  set. 

b.  oLars  with  north  declinations  greater  than  the  latitude 
will  exhibit  two  regions  of  near- zero  rate. 

G.  A  single  continuous  region  of  near- zero  azimuth  rate  occurs 
for  stars  of  declination  approximately  equal  to  North. 

d.  The  relative  duration  of  the  regions  of  near-zero  I'ate 
Increases  with 'declination  from  zero,  for  stars  with  de¬ 
clination  equal  to  latitude,  to  100^  fo'*  stars  with  de¬ 
clinations  approximately  equal  to  90*^  Noiwh. 

AZIMOTH  ANGULAK  ACCELERATION 

To  find  angular  acceleration  in  the  azimuth  coordinate,  differen¬ 
tiate  Eq.  (8)  with  respect  to  time  and  obtain,  after  simplification, 


u)  sm  20 


[ 


sin  ^-tanS 


+  [tdn Sc®s«^Oi*co^<ntj- 2 cos lotj  (lo) 

Equation  (lO)  set  equal  to  zero  may  be  solved  for  the  time  corres¬ 
ponding  to  the  raaximura  value  for  angular  >yelucily  in  azimuth.  The  values 
t  -  o,  G  =  0*^  or  l80°  are  found  to  be  mathematically  and  intuitively  cor¬ 
rect  .  Accordingly,  the  azimuthal  aiv/ular  velocities  previously  calculated 
for  t  =  o  are  the  m£ixlmum  rates  (as  Indicated  on  Fig.  2)  for  stars  pasoimig 
tangent  to  a.  10°  radius  circle  about  the  zenith. 
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Maximum  azimuthal  acceleration  is  found  in  the  vicinity  of  tut  = 

7.5°  either  side  of  star  transit.  By  finding  the  slope  Between  values 
of  0  calculated  for  cut  =  5°  =  10°  ,  a  value  e  =  4. 28- 10"° 

degrees  per  second  per  second  was  computed  for  the  maximum  accelera¬ 
tion.  The  velocity  at  this  time  is  6  =  1.55*10"’^  degi-ees  per  second. 

ELEVATION  ANGULAR  POSITION 

Using  Eq.  (30)  (reference  4,  page  35 )»  we  write  for  the  elevation 

angle 

sin  lr\  s  cos(^-l)»2cos  ^c6s£  (ii) 

sa  C  -  D  sin*  ^  (12) 

Case  1 

4>  =  390 

5=  490 

C  =  cos  ($  -I  )  =  cos  (39°  -49°) 

=  cos  (-10°)  =  0.98481 
D  =  2  cos  <i>  cos  S  2  cos  39°  <^os  49° 

=  2 '0,77715  *0.65606  -  1.01971 

Case  2 

$  =  39° 

S  =  290 

C  =  ('.ns  10°  =  0.98481 
D  =  2 -0.77715 *0.07462  =  1.35942 
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star  rise  and  star  set  occtir  et  zero  elevation.  From  Eq.  (12)^ 


this  occurs  at 


C  »  D  sm 


2  Hit 


.  2  (rtt  «  Cos(<»-8)  _  cos^co9B’¥s\n^s\^^_ 

I  ^  Zco&^ctsS 


ojt  =  2 


.  H  r  I  ++«*>  ♦  to«  4 iHt 

— j 


Case  1 


=:  2 


4- "tan  49^1  ^ 


s  X  sm 


.  -*»  rs4=  d.@09T8  »,  A  a£.s7ql^ 

•  AA  !■■*  I  I  MwiiiM  Miriii  riMi— *w— «aA  ST  w  “  -  . 


=  2  iw' 0.^8273  =  2(± 79*20')  ®  ±  158®-W 


CaoG  2 


-  2  s’m  ^  “  j 

=  2  0.724.A3]''^ 


s  2 


=  2  s5Yf'o.Mll+  =  2C±58*2o')»  i 


The  maKlraum 


elevation  sjigle  we,s  originally  chosen  as  80°  for  both 


k 


ELEVATION  ANGULAR  VELOGIUy  AKI)  ACCELERATION 

Talce  the  derivative  of  Eq.  (l2)  to  obtain  the  elevation  angular 
velocity  as  follows:  ■; 


COS  k  dh.  s 


«  D  u>  sm  ^  COB 

2.  2 


(16) 


Ji...  D  m  sm  ari: 

d-fc  *'  “ 


(17) 


^  4» 


p  u)  sin 


[4  -  (2  C  -  D  P  MS  *** 


(18) 


•  •• 

For  pTxrposes  of  numerical  evaluation  of  h,  h,  and  h  as  a  fimction 
of  time,  the  alternative  form  may  be  used  for  elevation  an.gular  rate 
as  given  by 


4k  «  P  u>  sin  ast 
"  2  co*K 


(19) 


Figure  h  shows  el.e’/utl'''n  angle  and  angular  rate  as  a  function  of 
time  for  both  Case  1  arid  Case  2. 


Elevation  angular  acceleration  is  foimd,  b;y  taking  the  derivative 
of  Eq.  (19).  Thus, 


ifH  =  h 


-  »  -T  .P  cos  h  CoS  b)t  -f*  sm  cot  s'm  It 

2eos*k 


(20) 


11 
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The  points  of  maximum  and  minimum  velocity  are  found  by  setting 
acceleration  equal  to  zero.  Accordingly, 


id  co«  K  cos  u)t  4*  sivi  toi:  siti  h  a  q 


(21) 


(d  00$  In  06$  = 


p  0)  siv\  V\ 

2  Qojs  1^"“ 


(22) 


COS  cot  ^  D  S\16  H 


sm^  w*t  Z  cos^  K 


(23) 


«V'  V  9  1 


- - L—  S  cos  4)  06S  S(C0S  to»t - 5^) 


cos  (ot  ‘ 


(24) 


Trial  of  values  in  the  heighborhood  of  ut  -  45*^  were  made  using 
Eq.  (24)  as  the  measure,  with  the  result  for  Case  2  that 


(dtj  ,  ^  i  4*7, is  degrees  (25) 

ImAX  n 

and,  substituting  into  Eq.  (12),  Eq.  (19),  and  Eq.  (20), 


t  50-23 

^  0.1943 


degrees 


degree.^  per  minute 


o 


(26) 

(27) 
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Tlae  Kinimum  value  for  angular  velocity  in  the  elevation  coordinate 
is  zero  at  transit  (ojt  ~  0,  h  =  80°). 

A  region  of  near-aero  elevation  r^ite  occurs  near  ma'iclmuni  elevation. 

f&xlraum  acceleration  (positive  rate  of  change  of  velocity)  in  the^ 
elevation  coordinate  occurs  at  the  horlaon  (h  =  O)  where  ojfc  =  i  158°40 
for  Case  1  and  cut  =  i  ll6°40'  for  Case  2.  Using  Eg.  (20)^  we  obtain 


ss,  —  u?‘  so$  ^  O0S  S  ees  (28) 

iViAiX 

aecel 

Case  1 

4>  =  39° 

8-  49° 

s=-l/l6  cos  39°  cos  49°  cos  158°40' 

lAeeel 

=  0.77715 •0.65606*0. 93144 
=  0.0297  degree  per  minute  pmr  minute 

Case  2 

4,  =  390 

5=  29° 

««  I 

^  I  moLjc 
I  AC«el 


=  -I/16  cos  39°  cos  29^  cos  ll6'^40' 

=  0.77715 -0.87462 -0.44070 
- - - 


a  V\ 


max, 

deed 


0.01906  degree  per  minute  per  minute 


PREQUMCy  DOMAIN  EQUIVALENT  OF  AZIMUTH  TIME  FUNCTION 

Figure  3  may  be  interpreted  as  a  time-domain  specification  of  the 
problem  of  azimuth  tracking  imposed  by  the  motion  of  the  hypothetical 
star  rejjresented  by  Case  2.  For  the  purposes  of  searro  system  design, 
the  frequency  domain  transformation  of  this  time  function  is  desired. 
This  transformation  is  obtained  by  Guillemln's  impulse  method  (1,5)? 
in  which  the  integrand  of  the  laplace  integral  is  converted  to  a  set 
of  impulses  in  order  that  the  integral  may  be  evaluated  without  numeri¬ 
cal  methods.  Steps  in  the  application  of  the  method  are  Indicated  on 
Fig,  5.  The  position  of  the  time  function  that  lias  between  star  rise 
and  star  set  is  replotted  in  Fig.  5a »  taking  the  new  origin  of  coor¬ 
dinates  at  the  time  and  azimuth  angle  of  star  rise. 

1 

To  obtain  the  set  of  impulses,  the  derivative  of  the  time  function 
is  first  plotted.  Fig.  5h,  and  approximated  by  a  broken  line  (this 
approximation  is  equivalent  to  fitting  the  t.lme  function  with  a  series 
of  parabolic  segments).  The  broken  line  approximation  Is  then  differ¬ 
entiated  twice.  Fig,  5c  and  to  obtain  a  set  of  impulses.  The 
change  in  units  from  radians/second^  to  radians/second^  associated  with 
the  talcing  of  the  third  derivative  is  offset  by  the  fact  that  the  in¬ 
dicated  strengths  are  related  to  areas  with  time  as  the  abscissa.  The 
impulses  shown  in  Fig.  5<3-  are  an  approximation  representing  the  Infor- 
matiqn  within  the  Interval  between  star  rise  and  star  set  ^vith  the 
assumption  that  a  continuous  function  is  involved.  That  is,  0,  6,  6, 
etc.,  have  the  same  value  immediately  prior  to  the  star  rise  point 
(t  =0)  as  immediately  after.  In  like  manner,  the  function  is  also 
continuous  at  the  end,  i.e.,  at  star  set. 

Api)endix  A  contains  the  details  of  the  matliematlcal  processing 
from  the  impuJ.ses  of  Fig.  5d  to  the  continuous  frequency  domain  equiva¬ 
lent  shown  in  Fig.  6. 

AZIMUTH  SERVO  SYSTEM  TRANSFER  FUNCTION 

In  the  Graham  design  procedure  (2),  use  is  made  of  a  steady- 
state  error  series  having  the  general  forifl 


where 

Kp  ---  position  error  constant 
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0(t) 


FIRST 

DERIVATIVE 


Fig.  5  -  Derivatives  of  the  time  function  in  the  approximation  process 
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RADIAN  SECONDS  (08  RELATIVE  TO  1  RADIAN  SECOND) 


Fig.  6  -  Azimuth  excitation  function  transformed 
into  the  frequency  domain 


=  velocity  error  constant 

=  acceleration  error  constant 

=  rate  of  change  of  acceleration  error  constant 

Anticipating  need  for  information  regarding^ the  input  motion  and  its 
derivatives,  expressions  for  Oi  ,  ©i  ,  and  ©i  were  derived  and  appear 
earlier  in  the  report.  Variations  in  design  predicated  upon  higher 
derivatives  than  the  second  (acceleration)  should  he  avoided  (6). 


Study  of  Pig.  6  leads  to  the  choice  of  the  servo  system  transfer 
characteristic  shown  in  Fig.  7-  The  low-frequency  asymptotic  slope 
of  -1  is  followed  hy  a  high-f?§quency  asymptotic  slope  of  -2  with 
unity  gain  at  the  break  frequency  (7)»  For  such  a  system,  the  steady- 
state  error  expression  is 


The  maximum  required  bandwidth,  i.e.,  maximum  value  for  ci:i3_  ,  is  found 
by  using  the  maximum  value  of  input  velocity  in  Eq.  (30)  •  For  an 
assumed  specification  of  l/U°  maximum  allowable  error,  a>i  =  0.084. 
Thus., ...a  bandwidth  of  0.1  radian  per  second  is  adequate  for  the  train 
axis.  If  the  tracking  accuracy  requirements  had  been  lower,  less  of  ai 
"overcoat  fit"  of  the  excitation  function  shown  in  Fig.  6  woiild  have 
resulted. 

AUTOMATIC  TRACKING  REQUIREMENTS  STARTING  FROM  REST 

Consider  a  radio  str  r  tracking  system  in  which,  the  antenna  feed 
point  is  nutated  to  porirh/.-  development  of  an  error  signal.  For  such 
a  system  it  would  be  possible  to  demand  of  the  servo  system  an  ability 
to  start  from  initial  conditions  in  which  an  error  in  0,  0,  0,  etc., 
exists. 

The  function  shown  in  Fig.  6  is  based  upon  zero  initial  errors. 

To  illustrate  the  effect  of  an  error,  consider  the  initial  conditions 
(where  the  minus  and  plus  subscripts  refer  to  the  instant  before  and 
the  instant  after  the  transient) 
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Fig.  7  “  Type  I  servo  system  open-loop 
asymptotic  characteristics 
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e_  =  e. 


as  in  Fig.  5a 


eoc.  =  u 


©+  )  6+  }  etc.,  same  as  in  Pig.  5a 

The  effect  of  Initial  conditions  is  to  introduce  an  Impulse  at 
t  =  0  in  Fig,  5g  which  becomes  a  doublet  (8)  in  Fig.  5'i*  The  strength 
is  established  by  the  step  at  t  =  0  in  Pig.  5b,  i.e.,  +0.473*  The 
approximation  of  the  transform  of  the  time  fiinction  is 

e*Gio)  s  Rt-  *id«  e<(.(A0  + 

i.f ro.4-73  ,  G sm2JF‘| ,  •  f  «  cos  2.fF  ll 
=-10*-{L-pr-+-pi— J+Jl.  p5  ■][  (: 

The  magnitude  of  Eq.  (3I)  is  given  by 


2  0.4736  sin  2.?  F  .  6^1^ 


The  zero  frequency  asymptotic  fvin.etion  is 


as  eoiapared  with  the  9,84*10^1  P  funcoion  given  in  Eq,  (All)  wiaere  no 
Gta,];ting  tran.r.ient  servo  requirements  are  involved.  If  such  a  servo 
system  is  desired,  the  excitation  runetion  represented  by  Eq.  (32)  may 
be  studied  in  the  sujne  mannc-c  as  Eq.  (a8).  The  transfer  characteristic 
should  be  a  Type  II,  that  is,  one  in  which  the  initial  asymptotic  slope 
is  -2,  followed  by  a  mid-frequency  asyiapuotic  sloije  of  -1,  and  ending 


6a75*_  [6.S7a  lo^'p 
L  w  J 


wit)i  a  high-frequency 


PP.0PO3ED  SERVO  SYSTEM 


slope  (9). 


Figure  8  shows,  in  block  diagram  form,  a  possible  servo  system 
arrangement  for  the  automatic  star  tra.cking  p'roblera  in  azimuth  in  which 
e.linri  nation  of  errors  '’ue  to  initial  transients  is  not  required  --  the 


SPACE  FEEDBACK 
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problem  and  system  of  Figs.  5,  6,  and  Y*  ‘I'be  error  expression  for 
the  system  shovm  is  substantially  identical  to  Sq.  (30). 

The  switch  shown  in  Fig.  8  provides  a  means  whereby  the  antenna 
may  be  held  in  a  fixed  position  until  tracking  is  desired.  The  intro¬ 
duction  of  a  fixed  voltage  corresponding  to  a  selectable  constant  an¬ 
gular  velocity  reduces  the  servo  tracking  problem  to  tracking  the 
difference  instead  of  the  total.  The  (£Q  of  the  equalizer  should  be 
approximately  10  oii  or  1  radian  per  second. 


CONCLUSIONS 

Star  tracking  with  an  elevation-azimuth  coordinate  system  has 
been  studied.  Time  functions  have  been  given  for  position,  velocity, 
and  aecelerat ion  components  in  both  axes.  Details  of  the  transforma¬ 
tion  to  the  frequency  domain  equivalent  and  the  design  of  a  suitable 
servo  system  for  the  azimuth  axis  are  included.  A  bandwidth  of  0.1 
radian  per  second  i's  shown  to  be  adequate  for  the  azimuth  servo  system. 
A  proposed  system  is  shown  in  block  diagi-am  form. 
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APPENDIX  A 


TIME  DOMAIN  TO  FREQUENCY  DOMAIN  TRANSFORMATION 

The  approximation  of  the  transform  of  the  time  function  shown  in 
Fig.  5  is 

4 

■»*  id*®  r  -j  <0  0.6*10 

@  (ji4  “  /;"r5  -0.213  e 

y**"®)  L  ^  4  4 

Jtc3.os-icr  -jtoa.SMO  -jo)  2.65*io 

+  l.8Se  4  9.08  e  -11.33  e 

.JI6255-I0'*  .  ...  ..+  - - .4 

- 11.33  e  4^.68  e 


•  41.15  e** 


«jtos*icr 

40.633e  -0,233  e  J  (Ai) 

Writing  the  exponentials  in  the  form  of  in-phase  and  quadratirre 
con^onents,  we  hav^g 

»  -^r-§'r^*233  +  0.633 (c&s  0.4*10^0) -j  Sm  0.4*10  to) 

4  1.85  (cos  a.05*(C^t0-j  siV2/)5*iO(o)49.08(cos2.3*10%--3  5'n2.5*l0 o) 
-  11.33  (cos  2.65'lO^cO-j  s<ii2,45'/6^C(j)-II.33(cos2.95’»/0^cO-^  Shi 2.95 

4 3.o8(cos  3.S  'Ilf (0  -j  si’n 3,l'/dfco)4l,S5(bo5  3.56'/olo-j  Sm3S^*16ftj) 

4  0.433(cos  S*lo^ (O-j  sm5'’lo^io)-a233^s5,4'llf 0  -j’5m  S;4*/ oj  ] 

mm 

(A2) 


Introducing  a  reference  angular  frequency  =  10”''*'  radian  per 
second  and  letting  F  ==  oo/co^  ,  we  have 


24 


where 


^  L/n  si'vi  /.«a  SiU  ^  Lts 

- ^  —  -...  -  .  .  -...  w.». 


-W 


+I.S?  siVi  2.05  F  +a.0«  sivi  2.5  F-11.33  2.45F HI. ?3  Sm  2.$5‘F 

=t=$.0?  SsV*.  3.1  r  'fl.lFSsV.  3.S5F-f  d.t33  Sm  SF-(3.Z33  SmS.tpj 

+  j  [^0.2.33  Cos  0-P  +O.433CdS0.tF  -i-l.SF  ddS  2.dS  F 

+  3.6Seos2.5F  -(1.33C6S  2.4>5‘F  -11.33  C&S  2.95  F 

•f9.0geos3.lF  ■Fl.g5c.os  3.55 F  +  0.633 eoS  5  F 

—0,233  Oos  S.4  F  (A3) 

4- 

= - ^ j  I -0.233  (siifi  O'F  +  slvi  S.toF) 

■hO.633  (sM  0.6F  +  Sm  S  F)  4  l.gsfsm  2.05  F  +  sm  3.55  F) 

+  9.0g(sm2.SFH*Sm  3.1  p)  -11.33 (^iv\2.65F-fSm  2,95F)1 

'Fj  ^--0.233  C<loS  O’F-f  CoS  S.^F^ 

+  0.633  (cos  0.4F+C0S5F)  +  l,g5(0©S  2.05 F  + CoS  3»55F 

+9,of  Ccd52.5F+Cos3.lF}-11.33(eos  «45F  •»•  cos  2.95F)1  | 

(A4)"^  J 

1 


0  Cjov^  —  s.vi  2,0  F '«■  J  css  2.S  ^  j' (— 

G  «  0.4G&  cos  2.gF  -J.266  cos  2.2  F 
-3.7  CO5  0.75F  --!g,l6  e^sO.3  F 
+  22.66  cos  0.15  F 


lA^; 


(A6) 


(si"  2-6  P  ■=®S2.S  F) 


The  magnitude  hecomes 


=  JO® 

F® 


Before  uiwnerical  avaliiation  is  made  o±'  £q.  (A7)  as  a  function  of 
frequency,  consider  its  nature  as  frequency  approaches  zeroi  The  expansion 


COS  X  =  I  -  -I-  ‘  • 

ms^y  tenuted  in  evaluating  G  for  F->0.  Tims, 

£i.-.|<3|=  &^UoAU(l-^^) 


=  J.2fe6  » 


”•  3.7  •  5.7^  —  |S,l&  }  (Aio) 

Substitution  of  Eq.  (Aiti)  Intr;  E'.!..  (a{)  yield;-' 


F->o 


4  •  IG'^  F 


^  2.74 -16^  2.S4 

F  ifii 


Equation  (All)  shows  that  the  excitation  function  frequency  domain 
equi valent ^  when  plotted  on  log-log  scales,  has  a  minus  one  slope  as 
the  frequency  approaches  zero.  This  minus  one  slope  emended  has  a 
value  of  Tonlty  at  a  frequency  F  «  m  »  (s  =  2.84*10^  or  cn  =  2.84. 

<ik)  10"^ 

The  series  expansion  of  the  cosine  functions  fails  to  facilitate 
similar  determination  of  the  high-frequency  asynptote.  Instead,  the 
cyclic  natiire  of  Eq.  (A7)  jmist  be  considered  and  a  determination  of 
the  peak  value  made.  The  asymptotes  are  tangent  near  the  peak  values. 

The  function  Itself  makes  cyclic  excursion  to  zero  or  minus  infinity 
on  the  log  plot.  Hence,  the  asymptotes  must  be  regarded  as  defining 
the  upper  limit  of  the  function. 

Repeating  Eq.  (A6), 

-G  sj-0.444>  eas  2.1  F  +  1.266  Cos  2.2  F  3.7  COS  <3.75 F 

+  \S.I6  COS  0.3  P  --22.66005  0.15^  (A12) 


may  be  -ften  as 

-<5  «  -0.466 CoS(^)lir+*I.266C6s(:|)2Tr  +  37<ios(^^)2-ir 


+  i1.i6  -  22.66  cos(^2Tr  (A13) 


where 


i-.aJL  .r-M  T  231  t-331  o«A  t  3^ 


using 


Considering  the  factor  2'tf  separately,  express  the  fractions 
the  least  common  denominator  as  follows; 


1  f  1  »  1  f  !•  __  >  1 

^  2.2'  0.75  0730  0.15 

become 


165 

752 


210  ,  616  ,  1540 
752  7^  "752 


3080 


(m4) 


(AI5) 


462 


The  nuafoer  of  cycles  of  each  frequency  In  the  camposite  period  is 
found  by  detejnaination  of  the  least  common  multiple  of  the  numerators 
of  the  ratios  (AI5).  Thus,  the  numbers 

165,  210,  616,  15^0,  3080  (AI6) 


lead  to 

56-165,  154-210,  15.616,  6.1540,  3.3080  (AI7) 

The  products  (AI7 )  show  that  the  number  of  cycles  of  each  cosine  term 
of  Eq.  (A7)  in  the  period  of  the  combination  is  56,  1|4,  15,  6,  and  3 
respectively.  The  relative  iagiortanee  of  the  terms  pi’ogresses  down¬ 
ward  as  frequency  increases.  Since  the  sign  of  the  two  low-frequency 
terms  are  opposite,  the  value  of  F  for  which  they  are  out  of  phase 
appears  to  be  a  region  where  a  maximum  occurs. 

In  the  general  case  great  difficulty  may  be  experienced  In  locating 
the  region  of  the  peak  value.  Setting  the  derivative  of  G  equal  to  zero 
results  in  a  sum  of  sine  terms  with  the  desired  value  of  F  Implicit.  A 
series  expansion  of  the  cosine  terms  in  G  prior  to  taking  the  derivative 
changes  the  problem  to  a  solution  of  a  high  order  polynomial  which,  since 
it  is  based  upon  a  limited  nmber  of  terms  of  the  series,  gives  only 
approximate  answers. 

A  simple,  easily  applied,  procedure  consists  of  setting  up  as 
many  rows  as  terms  (five  in  the  present  example)  and  pi'ovidlng  columns 
in  which  the  location  of  the  maximum  positive  and  the  maximum  negative 
value  may  be  indicated  for  all  cycles  of  each  cosine  term.  The  point 
of  favorable  coincidence  con  be  established  through  scanning  by  eye. 

This  was  done  for  the  exaiiq)le  using  the  following; 


Columns  per  cycle 

Term 

220 

cos  0.15  : 

110 

cos  0.30  : 

44 

COB  0.75  ■ 

15 

cos  2.2  ] 

165/14  ^11.78 

cos  2.8  ] 

28 


As  it  developed,  the  coaigjosite  function  showed  symmetry  about  the  h?tlf- 
cycle  point  euad  the  process  was  completed  using  330  columns.  Th.e  non¬ 
integral  columns  per  cycle  for  the  highest  frequency  term  caused  no 
parbicular  difficulty.  At  columns  110  and  550  a  value  36.^87  was  found. 
The  peak  occurs  at  column  330  with  the  value  37-92. 

Equation  (AC)  may  be  regarded  as  having  two  psirts.  The  denominator 
gives  rise  to  a  -3  slope  on  a  log- log  plot.  The  remaining  problem  is  to 
locate  this  -3  slope  so  that  tangency  T/lth  t'jie  cyclic  numerator  is  ob¬ 
tained.  The  numerator  wlH  have  a  -3  slope  (on  a  log- log  plot)  at  the 
points  of  tangency.  Pfoceed  by  setting 


i  (lo^S) 


(A18) 


■‘if' 

(10 Ci) 

Substitution  of  Eq.  (A6)  for  G  leads  to 


(AI9) 


-B  a  p ‘(0^^0.466*2.8  &m2.8F  +1.266*2.2 sivi  2.2. F 

+3.7  0.75  5mO.75F  +  l8.l6-0.3^wO.BF 

-  22,66  ‘0, 15  sm0.l5  FJ 

-+  IO^[o,4(»6  cos  2.8  F  -  ( .266  cos  2.2  F 

"B.7  cos  0.75F-j8,(6  cos  0.3  F 

+22.66  cos  0.i5Fj  (^o) 


a  (i.So48  sm  2‘8F  -2,7152  sm  2.2F  - 2.775 SiV»0.7SF 
^  -5>^48  sm  O.^F  +3.339  smO.sSP) +(0, 446  cos  2,8F 

W.U4  cos  2.2F-.3.7CC56.75FH8.J6C&S  0.3F+2244fi6Sa(5j^  (  A21  ) 
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The  equality  of  interest  is,  as  already  estahlished,  in  the  region 


p  =1.5^  =  20*0*  (A22) 

Trial  and  error  yields  the  approximation 

F  :=  63.5748  (A23) 

For  F  =  63.5748,  Eq.  (A8)  gives 

Fafo3.574«  ~  *  -I*  3.22. 

The  next  point  of  tangeney  occiirs  one  cycle  later  or  for 

F  =  63.5748  +  125.6637  =  189.2385 

at  which 

1  (A£6) 

The  high-frequency  asymptote  crosses  unity  gain  at 

f3  =  37.23745*10'^  =  .3723745*10^ 

F  =  71.94 

q)  =  7. 194 -10"  3  (A27) 

Therefore,  the  equation  of  the  higii-frequenny  asyjiiptote  is 

d^Qm)  „  37.2?7-II>  .  37. 2. 37 -10^ 

P  -*■<» 

The  low-frequency  asymptote  intersects  the  high-frequency  asymptote  at 


(A24) 

(A25) 


30 


r(7.l94  •(0 


-4 

3.62  -fo 


The  function,  as  given  by  Eq.  (A8),  and  the  asyB5)totes, 
Eq.  (All)  and  Eq.  (A28),  are  plotted  in  Pig.  6. 


as  given  by 


(A29) 


The  high-frequency  csymptotic  slope  of  -3  results  from  the  tliree 
derivatives  in  the  approx.lmation  process.  If,  instead  of  one  deriva¬ 
tive  before  the  broken- line  approximation,  two  derivatives  had  been 
taken,  then  the  high-frequency  asymptotic  slope  would  have  been  -4. 
Clearly,  since  the  final  high-frequency  slope  is  a  function  of  the 
total  number  of  derivatives  taken,  the  actxml  slope  lias  nd  significance. 
In  looking  at  the  entire  curve  of  Fig.  6  then,  one  must  keep  in  mind 
that  the  relative  accuracy  is  highest  at  the  low-frequency  end  and 
drops  to  zero  as  the  final  slope  is  reached.  Since  servo  systems  with 
greater  than  -2  slope  (through  tmity  gain)  become  unstable,  sufficient 
accuracy  for  the  present  purposes  is  provided  in  Pig.  6. 
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